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Abstract 

Up-going and Horizontal Tau Air-Showers, UpTaus and HorTaus, may trace Ultra 
High Energy Neutrino Tau Earth Skimming at the edge of the horizon. We show 
that such events even for minimal GZK neutrino fluxes could be detected by space 
telescopes such as the EUSO project: These Horizontal Tau Showers will track very 
long fan-like, multi-finger showers whose signature would be revealed by EUSO, 
OWL experiments. Moreover the additional imprint of their young secondaries (/i^ 
and 7 bundles with pair fiashes) might allow to disentangle their nature from the 
old UHECR secondaries in horizontal showers. Indeed at large zenith angles , the 
number of /x^ and secondary 7's from old, Ultra High Energy Cosmic Rays becomes 
comparable. On the contrary up- ward muon bundles from UpTaus and HorTaus may 
arise within a young shower with a larger gamma-muon ratio (~ 10^). Such a very 
characteristic imprint maybe observed by Crown detectors on Mountain, planes or 
balloons in space, as well as by 7 satellites in Space. We estimate the UpTaus and 
HorTaus rate from the Earth and the consequent event rate of /i^ bundles, whose 
fiux at 92°— 97° exceed the up-going muon fiux induced by the atmospheric neutrino. 

Key words: 

Cosmic Rays, Air-Showers, Neutrino, Tau, Muons 



1 Introduction 



The study of ultra high energy upward and horizontal r air showers produced 
by r neutrino interactions within the Earth crust has been considered in re- 
cent years as an alternative way to detect high energy neutrinos. The problem 
of r neutrinos crossing the Earth is indeed quite complicated because of the 
complex terrestrial neutrino opacity at different energies and angles of arrival. 
In addition, several factors have to be taken into account, such as the amount 
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of energy transferred in the - r lepton conversion, as well as the r energy 
losses and interaction lengths at different energies and materials. This makes 
the estimate of the links between the input neutrino - output r air shower very 
difficult. Such a prediction is further complicated by the existence of a long 
list of theoretical models for the incoming neutrino fluxes (GZK neutrinos, Z- 
burst model flux, E^'^ flat spectra, AGN neutrinos, topological defects). Many 
authors have investigated this signature, however the results are varied, of- 
ten in contradiction among themselves, and the expected rates may range over 
a few orders of magnitude (Fargion, Aiello, & Conversano 1999; Fargion 2002; 
Bertou et al. 2002; Feng et al. 2002; Bottai & Giurgola 2003, Tseng et al. 2003; 
Bugaev, Montaruli, & Sokalski 2004; Fargion et al. 2004; Jones et al. 2004; 
Yoshida et al. 2004). So far, the majority of the current studies on this topic is 
based on Monte-Carlo simulations assuming a particular model of the incom- 
ing neutrino flux. Some author have chosen ad hoc the crossing depth of the 
tau; other considered ad hoc maximal distance for the tau in flight. Most of 
the authors focus on the UpTaus tracks in underground detectors. In previous 
works we have presented a very simple analytical and numerical derivation 
(as well as its more sophisticated extensions) which takes into account, for 
any incoming angle, the main processes related to the neutrinos and r leptons 
propagation and the r energy losses within the Earth crust (see Fargion et al. 
2004a, Fargion et al. 2004b for details). Our numerical results are constrained 
by upper and lower bounds derived in more simple and tested approximations. 
We have shown how the effective volumes and masses are more severely re- 
duced at high energy and we included as a further constraint the role of the air 
dilution at high altitude, where r decay and the consequent air-shower may 
(or may not) take place. We showed (Fargion et al. 2004a) that our results 
give an estimate of the r air-shower event rates that exceeds earliest studies 
but they were comparable or even below more recent predictions (Yoshida et 
al. 2004). Secondly, we pointed out that the consequent /x"^, e^, 7 signature of 
HorTaus largely differs from that of horizontal UHECR backgrounds (Fargion 
et al. 2004b). 

In this paper we introduce in the calculation of the number of events an ad- 
ditional suppression factor related to the altitude at which air showers are 
observed. This guarantees the optimal extension and the largest flux for the 
shower to be detected at each observational h altitude. In particular we apply 
all our previous results to the calculation of the expected number of events at 
different altitudes in the atmosphere, with a particular attention to the EUSO 
telescope. EUSO is a space born detector that will be located at an altitude 
of about 400 km on the International Space Station (ISS), aimed at the the 
detection of UHECR. It consists of a wide angle UV telescope that will look 
downwards towards the Earth atmosphere to detect the fluorescence signals 
induced by UHECR showers in the atmosphere. Its aperture is such to cover 
a surface as large as 1.6 x 10^ km, therefore it will encompass AG ASA-HIRES 
and Auger areas. EUSO is scheduled to be launched in the 2009. Its ability to 
observe within the down ward atmosphere layer the long tracks of horizontal 
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showers may show their pecuhar opening, by geomagnetic fields, into three 
main fingers defined by their charges: (muon and electron pairs of opposite 
charge bent in opposite arcs and rectilinear 7 paths for tens and hundreds of 
km distances ). 



2 How to estimate the Ecirth-Skimming Volume-Mass 

To calculate the effective volume we assume that the neutrino traversing the 
Earth is transformed in a tau lepton at a depth x, after having travelled for a 
distance {D{6) — x). The column depth D{6) defined as / p{r)dl, the integral 
of the density p{r) of the Earth along the neutrino path at a given angle 
9 is shown in Fig. 2. The angle 9 is included between the neutrino arrival 
direction and the tangent plane to the earth at the observer location (9 — 0° 
corresponds to a beam of neutrinos tangential to the earth's surface) and it is 
complementary to the nadir angle at the same location. The probability for the 
neutrino with energy to survive until a distance {D{9) —x) is e-iD(9)-x)/L^ ^ 
while the probabihty for the tau to exit the Earth is e~^^'"'. On the other 
hand, as we will show in the next section, the probability for the outcoming 
r to emerge from the Earth keeping its primary energy E^-^ is e^^^^^i^^ (where 
Q-'-e/L^(i3) ^ Q-'-e/ir energy E^ > 3 x 10^'' eV). By the interaction length 
we mean the characteristic length for neutrino interaction; as we know its 
value may be associated to the inverse of the total cross-section axot = fee + 
(Tjvc, including both charged and neutral current interactions. It is possible 
to show that using the acc in the g-iD(e)-x)/L^ factor includes most of the u^- 
regeneration along the neutrino trajectory making simpler the mathematical 
approach (Fargion et al. 2004a). 
The effective volume per unit surface is given by 



where A is any arbitrary surface above the corresponding effective volume. 
For instance this expression has been first estimated for all the Earth. In this 
case A is just half of the terrestrial surface, due to the request of selecting only 
the upward direction. Under the assumption that the x depth is independent 
of and Ir, the above integral becomes: 
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where the energy of the neutrino Ei, has been expressed as a function of Et 
via the introduction of the parameter f] = Eu/E^^.. the fraction of energy 
transferred from the neutrino to the lepton. At energies greater than 10^^ eV, 
when all mechanisms of energy loss are neglected, r] = E^/E^-^ = E^/Et-. ~ 1.2, 
meaning that the 80% of the energy of the incoming neutrino is transferred 
to the newly born r after the v — N scattering (Gandhi 1996, 1998), (Fargion 
et. all 2004a)). Once the effective volume is found, we introduce an effective 
mass defined as 

= Pout^ (1) 



where pout is the density of the outer layer of the Earth crust: pout = 1-02 
(water) and 2.65 (rock). 

The expression of the effective volume in the most general case becomes (see 
Fargion et al 2004 for details) 

VrotiEr) ( L,(p){Er) \ f . „ „,„ 

] — - = -T^ — X / e i^i'CcIvEt)) sm9cos9d9 (2) 

A 1 ^r(0)(^T) ] J ^ ^ 



where the interaction length Ivr(/3) (shown in Fig. 2 and compared to Z^-) guar- 
antees a high energy outcoming r even if outcoming from a thinner Earth 
crust (see Fargion et al. 2004 for a more detailed discussion of Lt-^p)). 
We remind that the total neutrino cross section (Ji, consists of two main com- 
ponent, the charged current and neutral current terms, but the r production 
depends only on the dominant charged current whose role will appear later in 
the event rate number estimate. The interaction lengths L-rp, L^cc^ depends 
on the energy, but one should be careful on the energy meaning. Here we con- 
sider an incoming neutrino with energy E^., a prompt r with an energy Er^ 
at its birth place, and a final outgoing r escaping from the Earth with energy 
Et-^, after some energy losses inside the crust. The final r shower energy, which 
is the only observable quantity, is nearly corresponding to the latter value Et^ 
because of the negligible r energy losses in air. However we must be able to 
infer Er^ and the primary neutrino energy, E,j, to perform our calculation. The 
effective volume resulting from Eq. 2 calculated for a detector with a 1 km^ 
acceptance area is displayed in Fig. 3. 



3 From the GZK neutrino flux to Tau-Air Shower rate for EUSO 



After having introduced the effective volume we can estimate the outcoming 
event number rate for EUSO for any given neutrino fiux. The consequent event 
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rate for incoming neutrino fluxes may be easily derived by: 



where L^cc = ((7Jv,.r^)-^ = JEt^^ = ^ • ^^"^^ (to§^) " ^m-^ sec"^ sr-^ 
describes as a flat spectrum (a = 2) most of the GZK neutrino flux and as 
a linearly increasing spectrum (a = 1) the Z-burst model; is the density 
of the most external layer (either rock or water). The assumption on the flux 
may be changed at will and the event number will scale linearly according to 
the model. 

In Fig. 4 we show the expected number of event for EUSO where we have 
included the Earth's atmosphere and we have used L^(^p^, so that we may 
express the results as a function of the final r lepton energy. 
As one can see from Fig. 4, at energy E — 10^^ eV the general expected event 
rate is given by: 

/ ^uE, \ _ / Er 



Such number of events greatly exceed previous results by at least two orders 
of magnitude (Bottai et al. 2003) but are comparable, but below more recent 
estimates (Yoshida et al. 2004). 



4 The Tau air showers and secondary differential rates 



The differential number of events is given by 



^ q)^pi}ki (4) 

dQdt Li^cc 



We introduce now a differential expression of the number of events which 
allows to calculate the number of events as a function of the angle 6. We shall 
introduce a suppression factor that cares the finite lenght of the horizontal 
atmosphere (Fargion et all. 2004a). 

We can rewrite the expression of the effective volume given in Eq. 2 as a 
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differential volume for each arrival angle 9: 
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and we obtain the following expression for the differential rate of events 
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with r] = 1.2 and E^^ = eV. 

If we now integrate on the solid angle dfl (half side) we obtain the above 
formula multiplied by a factor 27r. 

Given the r number of events we can calculate the rates of /i, pairs and 7, 
which originates as secondary particles from the r decay. The number of muons 
is related to the total number of decaying pions and according to Matthews 
(2001) is given by 



Er 
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and we obtain finally 



{E, e)^N,- ^ ' 7 ^ ^ (9) 



ded(j)dtdA^ ' dEdQd4>dtdA 



We show in Fig. 5 the average differential rate of r's and the secondary , e 
and 7 bundles from the decay of r leptons. We find that the muon signal at the 
horizon, related to Earth skimming tau neutrinos is above 10~^^ - 10~^^ cm~^ 
s~^ sr~^. One should notice that the muonic background produced by atmo- 
spheric neutrinos (CR fi^ ^ Vatm ^ p"^) below the horizon approaches the 
value — 2 • lO^^'^ cm^^ s^^ sr^^, which is at least one order of magnitude 

lower than what we have obtained from our calculation for minimal GZK Vj. 
fluxes (see Fig. 5). 
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Moreover this muonic shower would have a significant 7 component, with a 
high number of photons - iV^/iV^ ~ 10^ (Cilhs & Sciutto 2001; see also Fig. 
15 in Cronin 2004) - because of the r decay channels into both charged and 
neutral pions. 

On the other hand, horizontal UHECRs will not represent a source of con- 
tamination for our signal because the horizontal 7's produced in the hadronic 
shower inside the atmosphere would be exponentially suppressed at large slant 
depth {Xjnax > 3000 g cm"^) and large zenith angles {d > 70°) (Cillis & Sci- 
utto 2001). Only /i^ can survive when propagating through the atmosphere 
at large zenith angles. Such muons would also be source of parasite 7 signal 
- due to the pair produced in the //^ decay in flight - but the gamma-to- 
muon- number ratio would be now approximately <^ 1. 

Therefore, this difference would allow to distinguish gamma-rich HorTaus 
from common horizonal gamma-poor UHECR events to reveal UHE earth- 
skimming i^T-'s . 

For a more precise approach to the calculation of the rate of events one has 
also to take into account that the number of events varies as a function of the 
height h at which the observer is located to detect the muonic, electronic and 
gamma shower. Therefore we introduce an additional factor 



where the parameter hg that we have introduced as 



hs = Rr (Er) sin e + sin 9 (10) 



defines the optimal height where the shower can reach its maximal extension 
at the corresponding energy Er- This is the sum of the height reached by the 
r in the atmosphere before its decay (i?^sin^^), where we have neglected r 
energy losses in the atmosphere, and the altitude reached by the secondary 
particles of the shower which is related to the parameter X^ax- Note that 
Rr = 4.9(^^/10^^ el/) km and X^ax/Pr = 5.7 + 0.46 ln(^^/10iW) km. Here 
we have considered the air density Pr — 1.25 x 10~^ g cm^'' constant and equal 
to the value at the sea level. At low energies (10^^ - 10^^ eV) the second term is 
dominant and pr can be considered as constant because the r lepton travels for 
less than 1 km before it decays. At higher energy (10^^ - 10^^ eV) the first term 
is dominant, and we can neglect the way the exact value of Pr changes with the 
altitude. In Fig. 6 we show the differential number flux (per unit area, energy 
and time) at different altitudes of p leptons originated by r decay. We have 
performed this calculation assuming an input GZK neutrino flux and that the 
Earth outer layer is made of rock (left panel) and water (right panel). We show 
how the expected number of events increase with the altitude of the observer. 
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In particular at h = 40-20 km, roughly the altitude of the atmosphere layer 
where HAS or HORTAUs are expanding to maximal power and where EUSO 
is following the UHECR showers, the differential number flux of muons and 
electromagnetic showers reach already its maximal asymptotic value. 
In Fig. 7 we display the differential number flux (per unit area, solid angle and 
time) of the secondary fi^ pair (from HORTAUs) as a function of the aperture 
6 which describes the line of sight below the horizon (i.e. 6 = 9 zenith — 90°). 
Again we have assumed an input GZK neutrino flux and an Earth outer 
layer made of rock. The two panels correspond to two different energies: 
= 10^'^ ev (Fig. 7, left panel) and E^ = 10^^ ev (Fig. 7, right panel). 
Again we show that at higher altitudes the muon, (as well as the gamma flux) 
from HORTAUs is high and that satellites orbiting around the Earth at a 
few hundreds of kilometers may better search Earth skimming tau-air showers 
secondaries, if they will turn their focus on the Horizons edges of the Earth, 
where they are mostly produced (See Fargion et all. 2004c) 



5 Conclusions 



Horizontal showers from normal hadrons (or gammas) are strongly depleted 
of their electromagnetic component because of the large slant depth {Xmax > 
3 X lO'^ g cm~^), while horizontal tau air showers are not. Indeed "young" Hor- 
Taus either of hadronic (67%) or electromagnetic (33%) nature at their peak 
shower activity are expected to have a large N^/N^ ratio, greater than 10^ 
(but with a characteristic energy ratio E^/E^j, < 10^). Old horizontal showers 
would have N^/N^ ~ 1. This difference would allow to distinguish and dis- 
entangle HorTaus from horizontal UHECR events, even in absence of good 
angular resolution, opening a new perspective in the UHE neutrino Astron- 
omy. The secondary fluxes of muons and gamma bundles made by incoming 
GZK neutrino fluxes and their HorTau showers, is well above the noise (by 
one-two orders of magnitude) made by up-going atmospheric muons The neu- 
trino signals at energies much above EeV may be even better probing the 
expected harder neutrino Z-Burst model spectra (0), (|23 l. 
The peak fluences we find in the /i and 7 component at the horizon (±5°) 
will give a signal well above the background produced by atmospheric z/'s. We 
have not discussed the albedo muons whose fluxes measured by Nemo-Decor 
experiments, 0^ ^ 10~^ cm~^ s~^ sr~^, are made mostly by single tracks (19). 
Because pair (or triple) bundle muons are much rarer ichn < lO'^^^u , ~ 
10^^^ cm~^ s^^ sr~^), the search and detection of muon bundles by GZK Hor- 
Taus at a minimal rate of 10~^^ cm~^ s~^ sr~^ (over an area of 10^ m^) will 
lead, in a year, to about 30 muons possibly clustered in five - ten multiple 
bundles. These events will be reinforced by hundreds or thousands of associ- 
ated coUinear gamma flashes. A detector with an area of few tens or hundreds 
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of square meters pointing to the horizon from the top of a mountain would 
be able to reveal the GZK v^- — t young showers (lori, Sergi & Fargion 2004). 
The characteristics of a prototype twin crown-like array detector to be placed 
on mountains, balloons, or satellites (Fargion2001a) will be discussed in detail 
(Fargion et. all. 2004c)). The simultaneous sharp 7 bundle at 0^ ~ lO^^^lO^^^ 
cm~^ s~^ sr~^ and the "burst" of electron pair at (pe+e- ~ 10~^ cm~^ si^^ 
would give evidence of unequivocal r signature. 

It should be reminded that the neutrino interaction enhancement by TeV new 
Physics would produce also an increase of hundreds or thousands time in Hor- 
Taus beyond a mountain Chain (like Auger) than standard weak interactions 
would do. Therefore Auger (a)must soon detect the Andes Shower Shadows 
toward the far away west side because of their absorption inside the moun- 
tain; (b) Auger must reveal the absence (or) the birth of young (Fargion, Aiello 
Conversano 1999), (Bertou et all. 2002), HorTaus created by tau born inside 
the the mountain themselves; they will be more abundant if New Physics at 
TeV is inducing larger neutrino-nucleon interactions. In a few years Auger 
might anyway be able to observe the expected GZK neutrinos inducing Hor- 
Taus at EeV energies. In conclusion we showed that an orbiting telescope such 
as EUSO experiment will be able to see at least half a dozen of events of Hor- 
Taus mainly enhanced along the Continental Shelves or Mountain edges. To 
conclude we want to remind that inclined-vertical PeVs r air showers (Up- 
Taus) would nearly always be source of 7 "burst" surviving the atmosphere 
opacity. These sharp UpTaus (with their companion HorTaus above and near 
EeVs) might be observed by satellites as brief Terrestrial Gamma Flashes 
(TGF). Indeed we identified a possible trace of such events in BATSE record 
(taken during the last decade) of 78 upgoing TGF possibly associated with 
galactic and extragalactic UHE neutrino sources (Fargion 2002). 
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Fig. 2. Left: Column depth as a function of the incoming angle having assumed the 
multi-layers structure given by the Earth Preliminary Model. Right: Comparison 
between -Z>r(/3) and Ir for rock and water. As one can see from the picture, -Z^r(/3) is 
shorter than Ij- at energies above 10^'' eV, thus it corresponds to a smaller effective 
volume where r's are produced while keeping most of the primary neutrino energy. 
The energy label on the x axis refers to the newly born tau for ^r(/3) • ^u-r 




\o'^ \f ^ ' \- ' ^ ' b ' ' * ' ^^ ' ' 

10 111 II) 10 111 10 10" 11) 

Filial Tau Air Shower lincrgy IcVj 

Fig. 3. Effective mass for UpTaus and HorTaus per km square unit area including 
the suppression factor due to the finite extension of the Earth's atmosphere (an 
horizontal length of 600 km). The curves obtained (the red line for an Earth outer 
layer made of water and the pink line for the rock) are compared with a simplified 
model of the Earth, considered as an homogeneous sphere of water (green line) and 
rock (blue line). Here we used the interaction length -^^t(/3) the volume is ex- 
pressed as a function of the final tau energy. Note that above 10~^ km^ the effective 
mass-volume in the energy range 3 x 10^^ - 10^^ eV is larger than the volume of 
the atmospheric layer, whose ability to convert downward neutrino in observable air 
shower is negligible. Only horizontal neutrino interaction in the atmosphere may be 
detected, at a much lower rate than the HorTaus ones. 
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Fig. 4. Number of EUSO Event for HorTaus in 3 years record as a function of the 
outgoing lepton tau (iT(/3) as interaction length), including the finite extension of 
the horizontal atmospheric layer. At energy Ej- = 10^^ eV, the event number is 
Nev = 3.0 {(t)yE^/bQ eV cm'^ s'^ sr'^) for the water and Nev = 6.0 {cfiyE^/bQ eV 
cm~^ s~^ sr~^) for the rock. The resulting number of events has been calculated for 
an initial GZK neutrino flux: oc E~'^ . 




Fig. 5. Left: The differential number of event rate of r leptons (HorTaus) for an input 
GZK neutrino flux. As in previous Figures we are assuming that r's are escaping 
from an Earth outer layer made of rock. Note the discontinuity at ~ 1 rad, due 
to the corresponding inner terrestrial higher density core (see Fig. 2). Right: The 
differential number of event rate of the secondary muons produced by the decay in 
flight of r leptons in the Earth's atmosphere: HorTaus. As in previous Figures we 
are assuming an input GZK neutrino flux and that r's are escaping from an Earth 
outer layer made of rock. 
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Fig. 6. Left: The differential number ffux (per unit area, energy and time) of the 
expected number of events of fj, leptons from r decay at different altitudes as a 
function of the energy of the incoming neutrino E;^. Here we have assumed an input 
GZK neutrino flux and an Earth outer layer made of rock. Right: The differential 
rate of the expected number of events of fi leptons from r decay at different altitudes 
as a function of the energy of the incoming neutrino Ei, . Here we have assumed an 
input GZK neutrino flux and an Earth outer layer made of water The higher the 
observatory the larger the flux, and for h = 400 km it can even exceed the "noise" 
due to the atmospheric upgoing muons {^f_i,atm ~ 2 x 10~^^ cm~^ s~^ sr~^). 




Fig. 7. The average differential number flux (per unit area, solid angle and time) of 
secondary muons from HorTaus at different altitudes as a function of the aperture 9 
which describes the line of sight below the horizon (i.e. 9 = 9zenith ~ 

90°). The flux 

has been calculated for two values of the energy of the r: Er = 10^^ eV {left panel) 
and Er = 10^* eV (right panel). We have assumed an input GZK-like neutrino 
flux and an Earth outer layer made of rock. Note the discontinuity of the angular 
spectrum in the left panel, due to the sharp density contrast of the Earth at ~ 1 
rad and the asymptotic behaviour for h ^ 1 km. Because the Earth is nearly 
opaque to EeV i^r, only Earth-skimming neutrinos nearly horizontally are visible. 
The discontinuity in the spectrum at ~ 1 rad occurs also at energies as large as 
10^*^ eV [right panel), but at a lower flux, thus it is not included in this figure. 
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